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Abstract 

We derive a family of matrix models which encode solutions to the Seiberg-Witten 
theory in 4 and 5 dimensions. Partition functions of these matrix models are equal to 
the corresponding Nekrasov partition functions, and their spectral curves are the Seiberg- 
Witten curves of the corresponding theories. In consequence of the geometric engineering, 
the 5-dimensional case provides a novel matrix model formulation of the topological string 
theory on a wide class of non-compact toric Calabi-Yau manifolds. This approach also 
unifies and generalizes other matrix models, such as the Eguchi-Yang matrix model, matrix 
models for bundles over P 1 , and Chern-Simons matrix models for lens spaces, which arise 
as various limits of our general result. 
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1 Introduction 



Finding the solution of Af = 2 supersymmetric gauge theories by Seiberg and Witten in 
terms of associated families of hyperelliptic Riemann surfaces, the Seiberg- Witten curves 
[HE], was a significant development in theoretical physics. The so-called Seiberg- Witten 
theory unifies many branches of physics and mathematics. From the physics perspective 
generalizations to theories with various gauge groups and matter contents were studied, 
as well as their relation to and the embedding in string theory. From the mathematical 
viewpoint the Seiberg- Witten solution gave important insights into the topology of four- 
manifolds. Also the solubility of Af = 2 theories turned out to be related to the underlying 
integrability, and the relations to various integrable systems were found. The literature on 
all these developments is immense, and the good starting point might be to consult the 
following reviews 0, SI El El [TJ [8] , as well as references therein. 

The holomorphic partition function Z is a central object in the theory. Its asymptotic 
expansion Z = exp ^] g=0 /l 29_2 jF 9 (a, A) in h is a generating function for gauge theory 
instanton numbers, which appear as coefficients of h, and inverse powers of the vevs a/ 
in the diagonal of the adjoint Higgs field. In aj 1 the JF 9 (a, A) have a finite radius of 
convergence as it has to be the case for physical terms in the effective action. In particular 
T = jF (a, A) is the prepotential, which determines the exact low energy gauge theory 
effective action up to two derivatives. The jF s (a, A) for g > multiply gravitational 
couplings of the form R 2 + F^~ 2 , where R + and F + are the self-dual parts of the curvature 
and the graviphoton field strength respectively. 

The fact that JF (a, A) can be calculated from periods of a family of hyperelliptic Rie- 
mann surface S(a, A) over a meromorphic one form differential dS was the main insight 
of [U [2]. However this was argued using global consistency conditions of the low energy 
effective action and not from the microscopic action itself. The latter argument was pro- 
vided by Nekrasov and Okounkov (9j dO]. In [9] Nekrasov developed a direct instanton 
calculcus and used localization techniques to determine Z to all orders in h as sums over 
two-dimensional partitions labeling the instanton configurations. In the thermodynamic 
limit the limiting shape of the partitions approaches the Seiberg- Witten curve E(a, A) and 
in [JO] it was proven that the periods of the latter over dS reproduce T to all orders in a;. 

The curve E(a, A) and the differential dS are the defining data of the Af = 2 supersym- 
metric theories in the sense that Z can be reconstructed from them. In particular one can 
view S(a, A) as the spectral curve of a putative matrix model and the A-cycle integrals of 
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dS as the fixed filling fractions and use the recursive solutions [TT] of the loop equation to 
reconstruct the T g . 

In this paper we find the explicit 1-matrix model for SU(n) Seiberg-Witten theory, 
whose spectral curve is the Seiberg-Witten curve. The strategy we use is to represent 
partitions hy N x N matrix integrals in a way recently proposed by Eynard |14| . which 
explicitly reproduces the Nekrasov partition function in the large N limit. The matrix 
model which we find reads 

Z M = [ OTe-^^W, VM = l[dx l l[(x t -x 3 ) 2 , (1) 

J Mat nN i iK j 

where M e Mat nN is nN x nN matrix and the measure VM involves the ordinary Van- 
dermonde determinant. In the large N limit the potential is given by 

n 

V M {x) =tx + 2^({x- ai) log(x - ai) - (x - a,)) , 
1=1 

where a; are vevs of the Higgs field, while t in the linear term encodes in particular the 
scale A. 

More precisely the fact that it is possible to write the Nekrasov partition function 
as a matrix integral is a consequence of its close relation to the Plancherel measure on 
partitions. The Nekrasov sum for SU (n) theory can be regarded as a generalization of the 
Plancherel measure to the case of n sets of partitions; equivalently it can be written as sum 
over one set of the so-called blended partitions. The properties of the Plancherel measure 
have been known for a long time; in particular it is known that in the thermodynamic limit 
it leads to the smooth limiting shape of large partitions known as the arc-sin law [H2 [T6] . 
This was used by pEU] to show that analogous limiting shapes arise for the case of the 
Nekrasov partition function, and that they encode the geometry of Seiberg-Witten curves 
of the corresponding J\f = 2 theories. Then, recently Eynard demonstrated how to rewrite 
the ordinary Plancherel measure, as well as its g-deformation, as a matrix integral [H]. In 
this paper we use similar methods to derive the above matrix models for SU (n) theory. 

We also derive a matrix model for a 5-dimensional theory compactified on a circle. We 
find that it it reads 

Z 5d = [ VMe~T^ v ^ M \ VM = Y[d Ul l[(2smh^^y, (2) 

JMat nN ■ {<j V 2 J 

which is related to a deformation of the 4-dimensional result. Now the measure is given 
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by the deformed Vandermonde determinant, while the potential reads 

n 

V 5d {u)=tu + -u 2 + 2j2^2{e u+ai ). (3) 

1=1 

As is well known, e.g. from the geometric engineering of gauge theories [18] , the partition 
function for 5-dimensional theory is equal to the partition function of topological strings 
on appropriate geometry. Therefore our matrix models also provide a new formulation of 
topological string theory, in the large radius limit of geometries which admit a limit to 
SU(n) theories. In particular, in the text we will provide matrix model expressions for 
various non-compact toric Calabi-Yau spaces. 

A very important feature of the matrix models which we derive is the fact that in various 
limits they reduce to other well-known matrix models. In particular, the 4-dimensional 
model (which itself is a limiting case of the 5-dimensional model) is a direct generalization 
of the Eguchi-Yang matrix model [T9l 120] and reduces to it for n — 1. The 5-dimensional 
model, also for n = 1, leads immediately to matrix models for line bundles over P 1 [2"T] . 
On the other hand, in the so called orbifold limit which involves t = and suppression 
of Li 2 terms, the 5-dimensional model becomes just the quadratic matrix model with the 
deformed potential, which is the Chern-Simons matrix model for lens spaces postulated in 
[22], [23] and analyzed in [SUES]. These other seemingly unrelated matrix models turn out 
to be just various corners of a single matrix model landscape. 

The paper is organized as follows. In section [2] we recall a necessary background on 
Seiberg-Witten theory and the Nekrasov partition function. In section[3]we derive a matrix 
model for 4-dimensional Af = 2 SU (n) theory and show that its spectral curve coincides 
with the Seiberg-Witten curve of SU(n) gauge theory. In section [H a matrix model for 
5-dimensional theory and its spectral curve are derived and analyzed. In section [5] we 
discuss its relation to topological string theory. In section [6] we recover other well-known 
matrix models as certain limits of our most general matrix model. Section [7] contains a 
discussion. 

2 Seiberg-Witten theory 

In this section we review the solution of the Seiberg-Witten theory in terms of two 
dimensional partitions, as well as the underlying family of curves, and set up a necessary 
notation. After introducing two-dimensional partitions and the Plancherel measure, we 
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discuss how its generalization leads to the Nekrasov-Okounkov partition function in 4 and 
then in 5 dimensions. 



2.1 Partitions and Plancherel measure 

A partition A of an integer number |A| is a set of non-negative, non-increasing integers 
A = (Ai, A 2 , . . .), A!>A 2 >...>0, ^A i = |A|. 

A partition can also be presented as a Young tableaux given by the set of aligned horizontal 
rows, each containing \ boxes. The number of rows of a partition is defined as the number 
of non-zero Aj. When considering statistical ensembles of partitions one can choose various 
measures. A very important one is the Plancherel measure defined as 

= f\( Xi - X i- i+j ) 2 = T[ I (4) 

J-l V j-i J J- J- hookiU) 2 ' K ' 

i<j J Dex v ' 

where hook(D) is the Hook-length of a given box in A, is the number of rows, and 
dimA is the dimension of the represenetation of the symmetric group corresponding to the 
Young- Tableaux. In the thermodynamic limit with respect to the Plancherel measure the 
partitions approach a limiting shape given by the so-called arc-sin law [15], [[6]. This shape 
is described by the arcsin function, hence its name. An important fact is that, after the 
rescaling leading to the smooth limiting shape, the limiting partition has a finite length. 
This phenomenon is called the arctic circle property. 

One can also consider a g-deformation of the Plancherel measure, denoted P q , which 
arises by replacing the integers in (pj by its g-deformation 

[h] = q~ h/2 ~ q h/2 , W = i[( q - l/2 -q i/2 ), 

so that 



i=l 



) ■ (5) 
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2.2 The Nekrasov partition function 

The partition function for Af = 2, SU(n) theory was derived by Nekrasov in [9] and 
analyzed further in great detail in [10]. It is a generalization of the Plancherel measure 
to the case of n sets of partitions. It therefore can be written down as a sum over a set 
of n partitions \f\ with I — 1, ... ,n labeling various partitions and index i > denoting 
i'th row of a given partition. This partition function depends on the scalar vevs a; via 
a-ik = ai — dk and it reads (9j CEO, [26], [27] 

z su(n) = A 2n ^Z x , (6) 

A=(A( 1 ),...,A(™)) 



z x = n 



aiu + h(j — i) 

where |A| = J2?=i |A (0 |. The sums in the above expression are performed over all possible 
partitions. To convert this partition function into a matrix model expression, below we 
consider the sums over partitions with at most N non-zero rows. In this case Nekrasov 
partition function is recovered in large N limit, where iV — ■> oo and 't Hooft coupling 
t = hN is fixed. More generally, we could introduce an independent number of rows iVj 
for each partition A®. However it would not affect the answer, because the final result is 
independent of N. 



2.3 Seiberg-Witten curve 

The solution of the Seiberg-Witten theory is encoded in the Seiberg-Witten curve. For 
SU(n) gauge theory this is a hyperelliptic curve of genus n — 1, which can be defined by 
the equation 

A n ( W + -)= P n (x), P n (x) =X n + ^x"" 1 + ...+«„, (7) 

for w, x G C, which is related to the familiar hyperelliptic form y 2 = P% — 4A 2n by the 
identification w = -^{y + P n ). The A- and B-periods of the differential 

d s-J- x ^L 

2m w 

of this curve correspond respectively to the vevs a>i and derivatives of the prepotential 

1 dT 
2m dai 
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Let us recall now how this curve emerges as the stationary point of functional [10] . 
which arise from the ADHM construction of the instanton moduli space [9]. The k in- 
stanton partition function can represented by contour integrals w.r.t. certain ADHM 
eigenvalues 4>i 

\e 1 + e 2 d4>! ] <f)u{(f)ij + e 1 + e 2 ) 



J2 Z * f II [ 2 7re ie2 P(0j)P(0/ + £l + e 2 )J JJ + ei )(0 7J + e ^ 



Here we wrote a more general expression with ei, e 2 equivariant parameters, which lead to 
the Seiberg-Witten solution under the identification h — e± — — e 2 (which we also assume 
in this paper). Moreover 4>u = <Pi ~ (ft J an d 

n 

P(x) = Y[{x-ai). (8) 

i=i 

One can now introduce the density of eigenvalues 

k 

p(x) = eie 2 ^5(x - 0/), 
i=i 

in terms of which the measure is approximated for small ei, e 2 by the saddle point method 

as 



A 2kn Z k ~ex V (— E A [p] 



^eie 2 

with p the saddle point of the functional 



E[p] = - J PV ' dxdy P ^l p( ^ 2 - 2 y dxp(x)\ogP(x). (9) 
Moreover, under the identification 

n 

p(x) = f(x) - ^2\x-ai\, 



the above functional is equivalent to the functional 



£[f] = \ f' V ' dxdyf"{x)f{y){x-yf{\og{x-y)-^ 



y<x 

P ' V ' dx dy (N + f{x)){N - f\y)) \og{y - x), (10) 

x<y 



whose stationary point corresponds to the partition which arises as a limit shape /* of 
Young diagrams in ([6]). This limit shape partition for SU (n) theory is given by the function 
f(x) of the form 

n 

f( x ) = ^2fi(x-ai). 

i=i 
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In fact, the stationary point /* of the functional S[f] is related to the function <p(x) 



f m (x) = ILe<p(x), (11) 

and this is <f(x) which explicitly encodes the Seiberg-Witten curve. More precisely, ip(x) = 
<&(x + iO), where $ is a certain conformal mapping described in detail in [10], such that 
the so-called bands of $ determine the cuts of the Seiberg-Witten curve. The above 
construction was extended in [28], where the limit shape equations where cast in the form 
of the eigenvalue distribution for all classical gauge groups. 



2.4 5-dimensional generalization 

So far we considered the 4-dimensional theory. It can be regarded as a limit of a 5- 
dimensional gauge theory compactified on a circle of circumference (3. For finite /3 the 
quantities in 5-dimensional theory, such as the partition function, correlation functions, or 
the spectral curve, are trigonometric analogues of the corresponding 4-dimensional quan- 
tities. In particular the partition function reads [U [TO] [26], [27] 

Zgo* = J2 (M) 2 "'%A, (12) 

A=(A( 1 ),...,A(")) 

sinhfto^Af-Af^-,) 
ft 11 sinh^f» ;i . + ?' - i) 

where |A| = YH=i \^ I- Here and in what follows we often use the (quantized) values 
Pi = ai/h, and p ik = p, - p k . 

This theory is again characterized by a Riemann surface and a meromorphic differential. 
The corresponding limit shape can be found from the variational problem for the following 
functional [10] 

£p[f] = r dxdy (N + f'(x))(N - f(y)) log {- sinh ^f^), (14) 

J x<y P * 

which is a direct generalization of ([101) . In section 14.21 we will show that the correspond- 
ing curve can be found in a complementary way in the matrix model formalism, using 
the Migdal-Muskhelishvili formula for the resolvent. The relation of 5-dimensional gauge 
theory to topological string on non-compact Calabi-Yau spaces will also be discussed in 
section [5l 
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3 Matrix model for 4-dimensional Seiberg-Witten the- 
ory 



In this section we derive matrix models for 4-dimensional Seiberg-Witten theory. Our 
strategy is to introduce an auxiliary dependence on a size N of two-dimensional partitions 
in the Nekrasov partition function, and replace the sum over these partitions by a sum over 
N x N matrices using the techniques presented in [14] . This leads to a matrix model with 
a very complicated potential, which includes such terms as logarithms of the T-function. 
However the result we are after is rederived in the iV — > oo limit, which turns into the 
ordinary large N limit in the matrix model formalism. In this limit we get a matrix model 
with a smooth, (relatively) nice and (relatively) familiar potential. We also discuss how 
the Seiberg-Witten curve arises as a spectral curve of this matrix model. 



3.1 The Nekrasov partition function as a matrix model 

Our aim is to rewrite the partition function of the Seiberg-Witten theory ([6]) as a matrix 
model. For simplicity, we consider first SU(2) case with n = 2. We denote A = A'- 1 -', 
jj = \( 2 \ and introduce 

hi = Xi — i + N + pi, hi = ^ - i + N + p 2 , (15) 

where we can assume p\ « p 2 , so that 

ki > k 2 > ■ ■ ■ > krsr > pi > hi > h 2 ■ ■ ■ > /i/v > Pi- (16) 

In terms of the variables hi, kj in the above range the SU(2) partition function can be 
rewritten as 

z su(2) ^ ^ \^yjh i+ ki) f -q hi ~ h i pj ki ~ k i pj hi ~ h. y . 

The products in the above expression can be reorganized so that only non-trivial hi, ki 
with i = 1, . . . , N occur explicitly. Furthermore the form of ( fi~5l) leads for i > N to many 
cancellations. Therefore we get for the first product above 

u h -0 - ( n ^-^(fifi'^j^ 

i<j J l<i<j<N 1=1 j=l J 

N 



( n (^-M)(n7^ 



l<i<j<N 
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\(hi-pi)\ 



and the analogous expression for the second product. The third product involves differences 
(hi — kj) where the index i can be equal, smaller or greater than j. Similar cancellations 
as above lead to 

ft^ = (n— )( ft ^)=4nc*.-*.) 

l=J=l 1=1 2=iV+l 1=1 

n^h - w( n (fc-^Xn^) 

»<j l<i<j<N i=l v * ^ y 

n?^r - <-», + M n (^-^(n^) 

«>i l<j<i<N j=l v 3 

(17) 

respectively. 

Collecting all above contributions we observe that each difference (hi — hj), (fa — fa) and 
(hi — kj) appears twice. There are also terms with factorials involving differences with p\ 
and p 2 , which appear for both hi and kj on the same footing. We can therefore introduce 
one set of variables 

h=i,...,2N = (fa, ...,k N ,hi,..., h N ). (18) 
In terms of U the partition function can be written in a compact form as 

2N . 2nl ■ 

z»« ~ e n b-hm ll .-J*<;.-ri„ = 

h>l2>...>l N >P2>lN+l>->l2N>Pll<i<j<^N i=l K 1 V % 

E IT ('■- h)'U ,,. _„,)„(/. _-,)» ■ < 1£ » 

Zl,...,ZjV>P2>ZjV+l,-,i2JV>Pl l<i<j<2N i=l Kl l ± ) ' t>A)- 

which is normalized with respect to the constant factors like p\ 2 and N\ 2 ; they are irrele- 
vant for the further discussion and drop out anyway in correlation functions, thus can be 
safely skipped. In case when negatives quantities arise under the factorials in the above 
product, we replace these factorials by T functions, which provides the relevant analytic 
continuation. 

The expression (fT9l) already has a flavour of the matrix model, with the first product 
being the Vandermonde determinant. To write it explicitly as a matrix integral, we repeat 
the trick from [14] and introduce the function 



f(x) = — xT(— x)T(x)e 



„ p -VKX 

-iwx " c 



sm{7tx) 
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which has simple poles at all integer values of the argument. Upon integration along the 
contour C encircling [pi,oo[ part of the real axis, this function can be used to pick up 
all integer values of hi G [pi,oo[. Similarly, one would have to introduce another contour 
C encircling [p 2 , oo[ half-line to pick up relevant values of fcj. In fact, we can use the 
same contour C for both hi and kj, because additional kj &)pi,P2[ picked up by the latter 
integral contribute zero (for such kj, T(kj — P2 + 1)~ 2 has a double zero which cancels the 
simple pole of /(%))■ This leads to the following expression for (fT9l) 



» 2N 

z su ^ =i dxi ... dX2N n (xi - x 3 f n 



f(xi)A 



2nx, 



1 T(xi -pi + l) 2 T(xi -P2 + l) 2 

l<«<j<2iV i=l K F ' K ^ 1 

PMe-^W, (20) 



lMat 2N (C) 

where M G Mat 2 N(C) is 2N x 2N matrix with eigenvalues in a set C. Choosing real 
values of pi we can in fact reduce the contour to the real axis in the large iV limit. It is 
convenient to rescale x 1— > x/H, and in terms of this rescaled variable the potential reads 

-V(x) = -(2nlogA)ar-log/(|) + ^ logT^r 1 - W + l) 2 = 

1=1,2 

= -(2ni og A)x - log (|) - log r( - 1) - log r(~) + ^ + 

+2 ( log (xr 1 - Pl ) + log r (xh- 1 - p,)) . (21) 
1 

This expression looks quite complicated. Nonetheless, we will see that in the large iV limit 
it simplifies enormousely. 

Now it also becomes clear that in generalization to arbitrary SU{n) theory we need 
to introduce U for i = l,...,nN, which is a concatenation of strings of TV eigenvalues 
corresponding to matrices corresponding to each partition A^. In the last sum of (I2TT) we 
also need to take the range of I from 1 to n. 

In what follows also the derivative of this potential plays an important role. It involves 
the logarithmic derivative of the T function, i.e. the digamma function 

m = ± xlog r {x) = m. ,22) 

This function shares a number of interesting properties, for example 

00 00 „ 

ip(x + 1) = -7 - V (— ) = log X + — - V 



'x + n n' 2x ^— ' 2ix 2t 
n=i i=i 
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where 7 is Euler-Mascheroni constant and B 2 i are Bernoulli numbers. From this we also 
deduce 

00 R 



lf)(x) = --^ +log x- 
From these expressions we find 



2x ^ 2ix 2i 



1 R 

V\x) = -2nniogA + ^(-- = i-— + 21og(|-^)-E 



xh 1 — pi h ^—^i{xh 1 — pi) 

h . .A fi 2i B 



-2n(log ft + ft log A) + £ (— + 2 k*(* - «,) - £ jg^) . (28) 



It is important to stress that the potential flzlj) . and so of course its derivative above, 
explicitly do not depend on the size of matrices N. For finite N, similarly as discussed 
in [H], there could be at most exponentially small corrections in N, which vanish in the 
large N limit we consider in what follows. 

The above expression can also be integrated and yields 

V{x) =tx + 2j2((x- a,) log(> - a,) - {x - a,)) + 0(H), (24) 

where coefficients of linear terms are encoded in the constant t. We stress that so far we 
considered a finite value of N corresponding to the number of rows in partitions which 
appear in the Nekrasov partition function (jOj). Therefore the original Nekrasov partition 
function is automatically obtained in the large N limit for the matrix model expression 



N -»• oo, h -> 0, HN = const. (25) 

In this limit all the terms of order 0{K) vanish, so we conclude that the SU(n) Seiberg- 
Witten theory can be formulated in terms of the 1-matrix model with the potential 

n 

V M (x) =tx + 2^~2^(x- a ; ) log(x - ai) - (x - a t )\ , (26) 
l=i 

where t in the linear term encodes in particular the scale A. This potential is shown in 
figure [H For SU{n) theory it clearly has n minima to which we can associate n cuts, giving 
rise to the Seiberg-Witten spectral curve of genus (n — 1). 



1 Some complication due to explicit N dependence in eq. (2.23) in [O] arouse due to higher tk couplings, 
which we do not consider here. 
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Figure 1: Matrix model potential for SU(3) Seiberg-Witten theory. Left: the 
potential corresponding to the exact expression (f2TT) for finite N; complicated 
spikes arise from logT terms. Right: the potential in the large N limit, given by 
(1261 ). Three minima give rise to three cuts and genus 2 Seiberg-Witten curve. 



3.2 Seiberg-Witten curve as a spectral curve 

Let us recall some basic facts concerning matrix models. For the matrix integral 

Z= [ Me-^W 



over N x N matrices, the eigenvalues are subject to the effective potential 

1 - 

S eff (X) = 2 1^ lAi - A, I + v (^)- 



h 

i<j i=l 



Replacing discrete eigenvalues by a continuous eigenvalue density r(x) 

N 

X 



i N r 

— h(\i) I— > / h{x)r{x) dx 
i=i J 



leads to the effective functional 

dx dyr(x)r(y) log \x — y\ — I dxr(x)V(x), (27) 



Seff[r] = [ 



and the spectral curve of the matrix model is encoded in the discontinuities of the resolvent 
which leads to the minimization of this functional. 
Let us now introduce R(x) such that 

R'(x)=r{x), (28) 
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and integrate (1271) by parts to get 

3 eff [R]= [ dxdy R ^ R ^ + [ dxR{x)V'{x). (29) 

Jx^y i x — y) J 

We focus now on the matrix model for Seiberg-Witten theory (l20l) , with a potential whose 
derivative is given in (l23l ). Let us stress that the matrix model (1201 still depends on both 
N (size of a matrix) and h. Upon extracting a spectral curve corresponding to the Seiberg- 
Witten theory we should again take the limit (1251) in which all terms in (l23|) proportional 
to h are neglected. The only terms which survive are 

V M '(x) = 2^1og(a;-a,) =21ogP(x), (30) 
i 

where P(x) = YYi=i( x ~ a i) coincides with ([8]). Plugging this into ( 1271) we get 

S eff [R}= [ dxdy R ^ R ^ + 2 [ dxR{x) log P(x). (31) 
Jx^y ( X — y) J 

We observe that this functional is identical to j9]) if we identify R(x) = p(x). From the 
analysis in [10] we know that the stationary point of (J9j) leads to the Seiberg-Witten curve 
encoded in the derivative of the limiting profile (fTTl) . In our case we also need to take the 
derivative (1281 ) of R(x) to get r(x) corresponding to the eigenvalue density. We conclude 
that the stationary point of (1271) - encoding the spectral curve of our matrix model - 
corresponds to the same Seiberg-Witten curve as the one derived in [10J. 

The above argument relies on the form of the functional ([9]) considered by Nekrasov- 
Okounkov. There is however a way to obtain the spectral curve directly within the matrix 
model formalism. In more generality this works also for the 5-dimensional matrix model 
which we introduce in the next section. One can therefore derive the spectral curve for 
the 5-dimensional theory according to the procedure presented in the subsection 14.21 and 
then recover the curve discussed above in the 4-dimensional limit. 



4 Matrix model for 5-dimensional gauge theory 

In this section we generalize the construction of matrix models for 4-dimensional gauge 
theories in section [3] to a canonical class of 5-dimensional gauge theories. This construction 
will be further slightly generalized in section [5l 
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4.1 5-dimensional partition function as a matrix model 

We now wish to introduce a matrix model for the 5-dimensional SU(n) theory with a 
partition function given in (fl~2l) . Similarly as before, to start with we consider SU(2) case 
with n — 2. We denote A = A* 1 ) and fx = X^ 2 \ which are partitions with no more than N 
rows, and similarly as in (fl~5l) introduce 

hi = Xi-i + N + p u h = + N + p 2 . (32) 

Moreover we assume pi « p 2 , so that 

k\ > k 2 > . . . > /c/v > P2 > hi > /12 • • • > hjf > p\. (33) 

We also introduce q = e~ 9s = e _/3ft and the standard g-deformed notation 

[h] = q- h/2 ~ q h/2 , [h]\ = f[(q- i/2 -q^). 

i=l 

The SU (2) partition function now takes the form 

with (hi), (kj) satisfying the condition ( 1331) . These terms can be rewritten analogously as 
in section 13 . 1 L For example 

nn 1 ^ =( n ^(n^), 

i<j i<j 1 Jl l<i<j<N i=l L 1 J 

and similarly for other terms. Therefore, introducing one set of variables 



h=i,...,2N — (ki, • • • , k N , hi, ... , h N ), (34) 

we get 



2N 



^ w \^ t MJ'-^w^B^> ,35) 

To write this expression as a matrix model integral we again take advantage of a few 
facts discussed in [H]. First of all, we introduce the following notation for the quantum 
dilogarithm 

oo 

g(x) = l[(l-x- 1 q i ). 

i=l 
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It vanishes g(q h ) = for h a positive integer, and at such points its derivative is 

g(1) 2 efa h g -Mfc-D/2 

^ W J "' Q h (l - ? fc )<7(<T fc ) ' 
Therefore the following function has simple poles with residue 1 for x = q h with h G N 



/(*) = 



2 _i2Li OKX _I!^L 



(1 - x)y/xg(x)g(x 
Let us also note that for x — q l 



= q 



-l(l+i)/ A _ 9(1) 



(jT^i - ^^-p(^->o g (^),og(^«)). 
In what follows we use the notation 

Xi = q k = e~ 9sk = e Ui , (36) 

so that Obi IS *^ 

cylindrical coordinate. Now we can write 



2N 



m\ynli q (l-2N)k 



z °r ~ e 11 v-*)'n lv ._ lhmi ,, )2 

(I,),,, „l<i<J<'2« i-l [lh PU'Ml'. P21-; 



1 ^ / ^-^-O AT 



(37) 



l<i<j<2jV 

X II — r? *-i U[ " 6 21 ° g9 



i=l Xs 



5(1) 

The deformed Vandermonde determinant can also be written as 

ed-AOE.iog*, Y[( Xi - x 3 f = J] (2sinh ig|gLZjg!gl ) 2 = A J. (38) 

j<j i<j 

With the deformed measure 

VM = A 2 q Yl dui = A 2 q e~ ^ JJ ^ 

i i 

we can rewrite the above partition function as a matrix model expression, which in a 
straightforward way generalizes to arbitrary SU (n) gauge group 



Z SU(n) = f pMe __L Tr y 5d (M) ; 
JM„ N (C n ) 



'M nN (C q ) 
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with the deformed integration contour C q is a circle of radius r e]l, q 1 [. The matrix model 
potential takes the form 

■2n log(/3A) 



-V 5d (x) 
9s 



9s 



+ (n — 1)N — ij logx — 2 logg(l) + logg(x) + log g(x x ) + 



1 Z7T I 

+ - logS + log(l — X) H logs (logx) 2 — 27T + 

2 5- s 2# s 

+ ]T [2 logp(l) - 2 log<7(arY0 - ( lQ g^ P ')) _ 1 log^g-ft^ 



2 logg 



(39) 



Similarly as in 4-dimensional case, this expression also simplifies in the large N limit. 
To turn the potential into a nicer form, we again differentiate it and then integrate. First 
we note [14] that 



log 0(2;) = VLi 

Qs — 

ys m=0 



where B m are Bernoulli numbers. The polylogarithm is defined as 

00 j 

Li m (x) = ^|-, 

and it has the following properties 



i=i 



Li'„(x) 
Li m (x) 
Li (x) 
Lii(x) 



Li n _i(x) 



-lr+Wx- 1 : 



for m < 



x 



1 — x 

— loefl — x). 



(40) 



The following useful relations follow 



g'(x) 
g(x) 



d 



x—(\ogg(x) + \ogg(x x )) 



1 t • / -1\ "msi 

> Lii_ m (x ) r 

ys m=0 

log x 1 x + 1 m 



g s 2 x - 1 g s 

Using the above facts, after some algebra one therefore finds the derivative of the 
potential 

T/ , f v _ 2nlog(/?A)+^((n-l)iV-i) 

^dW — r 



X 



x — qPi x 



2 1 1 2 - — ^ 

- log {(xq- pi Y^ - (xq~ Pl )^) + - ^h-2m(xq 



_ pl ,B 2m g. 



2m 



m=l 



(2m)! 



(41) 
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Let us note that log ((xg -P! )~2 — (xq^ 1 )^ is the q-deformation of the 4-dimensional expres- 
sion log (a: — ai) in (1231) . while other terms (apart from the linear in 1/x) are proportional 
to higher powers of g s . The expression which we obtained can now be integrated. In 
particular J ^log(x~ 1//2 — x 1 / 2 ) = — ^(logx) 2 — Li 2 (x), so that 

V 5d = t\ogx + ^(logx) 2 + 2^Li 2 (xe a; ) + 0(g s ), (42) 

i 

where the coefficient t encodes in particular the scale (3h and terms linear in a\. In the 
large N limit 

N — > oo, g s — > 0, g s N = const, 

we can again neglect all terms 0(g s ). We also note that introduction of the single parameter 
t is justified, as it involves both log(/3A) and g s N = const related to each other in large iV 
dualities precisely in this way. Substituting now x = e u , finally we get the matrix model 

z su(n) = f VMe -^v^M) = f ^ du .JJ (2sinh^^)V^^ y6d ^ (43) 

i i<j 

with the potential 

n 

V" ,d {u)=tu+ r ^u 2 + 2^U 2 {e u+ai ). (44) 

i=i 

The coefficient n in the quadratic term can be interpreted as redefining the coupling in 
front of the potential as 

9s 

9s ^ 9s = — 

n 

(and rescaling other terms in the potential by n). It is consistent with the same rescaling 
found in a related matrix model in [23], which we discuss in sectional The potential V 5d 
is shown in figure El 

4.2 5-dimensional spectral curve 

Finding the spectral curve for 5-dimensional theory also requires the 't Hooft limit, in 
which all terms in the potential (|4*T1) proportional to positive powers of g s vanish, so that 

n 

xVi d {x) — > -2Y,teg({xq-™)- 1 *-{xq-^)=-2\ogP q {y), 
i=i 

n 

for P q (y) = JJ[y-pi], 
i=i 
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Figure 2: Matrix model potential for 5-dimensional SU(3) Seiberg-Witten the- 
ory. Left: the potential in the large N limit, given by (1441 ). Three minima 
give rise to three cuts and the genus 2 Seiberg-Witten curve. Right: zooming 
around the rightmost minimum. 



with x = q v . This is indeed the trigonometric counterpart of the 4-dimensional case 
(l30l) . Accordingly the 5-dimensional curve will be the trigonometric counterpart of the 4- 
dimensional Seiberg-Witten curve, in agreement with the relation between 4-dimensional 
and the 5-dimensional variational problems (fTUl) and ([H]). 

Apart from arguments which rely on the form of effective functionals (which were also 
employed in section [3721) . we can in fact derive the form of the 5-dimensional spectral curve 
using purely matrix model technology, i.e. determining the resolvent and the eigenvalue 
density in the large N limit. This method could also have been be applied explicitly in 



section 13.21 however it is more convenient to consider the general 5-dimensional case and 
then get the 4-dimensional results in the gauge theory limit. For a general potential in 
the multi-cut case the resolvent uj(z) is given by the Migdal-Muskhelishvili formula [29) . 
It has the following form in terms of the periodic variable x 



where C = C\ U . . . U C n is a union of n cuts supported on intervals Ck = [x2k-i, x 2k\- 
The 2n endpoints of the cuts Xk must be suitably chosen. A set of (n + 1) conditions for 
these endpoints stems from the limiting behaviour of the resolvent at infinity lu(z) ~ 1/z. 
Another (n — 1) conditions can be chosen in a way which is appropriate for a given physical 
problem. In our case it is natural to fix (n — 1) filling fractions, which are determined by 




(45) 
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the A-periods around the cuts. 

To compute the resolvent (l4*5l) one can follow the solution presented in [21], where the 
case corresponding to n — 1 (which we discuss also in section 16.21) was analyzed in detail. 
This solution relies on the chiral ansatz and the arctic circle property of two-dimensional 
partitions which we discussed above. We assume that the eigenvalue density p(z) is non- 
trivial only in certain intervals, independently for each set of partitions summed over in 
the Nekrasov partition function. The parameters ai in the potential (1441) can be absorbed 
into the positions of the endpoints of these integrals, so that the resolvent becomes 



LU(Z) 



2ni 



n 



dx- 



c 



logx + t - 2 YJl=i l°g(l ~ xe 
x(z — x) 



nz - e-^fc- 1 
\ x - e-y™ 

\ k=l 



(46) 



To perform this integral one can modify the integration contour, so that it encircles the 
branch cuts of all logarithms and the pole at x = z. 



-y, 




Figure 3: Integration contour used to compute the resolvent for 5-dimensional 
SU (2) Seiberg-Witten theory. 



As an example let us consider how to derive the spectral curve for the SU (2) theory 
from our matrix model. We denote the endpoints of the two cuts as [2/1 , 2/2] and [1/3,1/4]. 
The integration contour in this case is shown in figure [3l The resolvent becomes a sum of 
the residue at x = z and the contour integral along the (infinitesimal) circle around x = 
(a part of the contour around the branch cut of log x) , as well as contributions from the 
integrals along parts of real axis from branch cuts of all logarithms. The latter are given 
by integrals of the form 

/ /I («> 

v(z -v)\/Yl i=1 (v -e-y>) 

They can be expressed in terms of elliptic functions of the first and the third kind, as 
explained in the appendix. In what follows we also use the notation introduced in the 
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appendix. After some simplifications, for SU(2) theory the resolvent (Pf6l) turns out to 
involve three integrals of the form (11711 evaluated at the points v = 0, oo, e~ V2 (the integral 
evaluated at e~ m vanish). For these points the arguments of the elliptic functions F and 
II are respectively 

e 2/3 (e 2/4 — e 2/2 J e 1/4 — e 2/2 

as well as 5 and A; 2 given in (1551) with a, b, c, d identified with e~ Vi . Using the addition 
formula ( 16T1 ) twice we can replace three elliptic functions F(y(v), k) by a single one with the 
argument 2/0,00,2/2 determined by using ( 1641) twice. Similarly, using ( 1621) twice we get a single 
function Il(?/o,oo,j/2 5 ^ but then in addition two terms (artanh 0o,oo) and (artanh Oo,oo,j/ 2 ); 
where Oo,oo and Qo,oo, y2 are determined from the formula (l63l) used twice. Altogether, the 
resolvent becomes 

t + 21ogz- 2 log(l-^ 2 ) -21og(l-ze w ) x/ILO 2 = e_K ) »i+«a+w+»4 
wte) = — — + v ' -te 2 + 

2 2 



+ 



VIL(*-e- 



+ e _, 3e _, 4 (n(y(0), 5(0), fc) + n(y(oo), 5(0), k) - n( 2/ (e^ 2 ), 5(0), fc) J + 
_i/4 — p ~ y3 / artanh 0o,oo + artanh ©0,00,2/2) \ 



n(y <5(z),A;) + 



(z-e-^J^-e-w) V v» u '°°'W> w> v /5( 2 )-i( ( j( z ) - l)(5(z) - fc 2 ) / 

We discuss now the behavior of this resolvent at infinity z — > 00, where it should decay 
like ~ I/2. Firstly, we note that the last term in the first line, the last term ~ e m F(yo t00iy2 ) 
in the second line, and the terms in the third line are of order ~ z. Therefore a sum of 
coefficients of these terms must be zero, which gives one equation for the endpoints of 
the cuts. Due to simple the same dependence on z, these terms in fact cancel altogether 
and drop out from the expression. Then, the only term which is constant at infinity is 
the first term ~ F(y 0oo ^ 2 )/(e~ y3 — z) in the second line. It vanishes if and only if the 
elliptic function ^(2/0,00,2/2) vanishes, which is possible if and only if simp = yo,oo, y2 — 0. 
One one hand, this provides another condition on the endpoints of the cuts, due to (1641) . 
On the other hand, then the terms in the fourth line II (7/0,00,3/2 = 0,5(z),/c), as well as 
(artanh Oo i0 o,j/ 2 ) (whose argument is proportional to sin cp = 1/0,00,2/2 Via (|63l) ) also vanish. 
Therefore the only term which could contribute to the eigenvalue density (or equivalently 
the spectral curve) is (artanh ©0,00) in the fourth line. Quite non-trivial cancellations lead 
to its overall prefactor being just 2/z which in particular does not depend on any e~ m . 
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Moreover, this term indeed has the square root branching as hoped for, which arises via 
the dependence on z through the argument of 6 0jO o 



^5(z)(l-5(z))(S(z)-k 2 )~ ^U Z 6 



-in • 



e-y* - z 



2 



Indeed some algebra reveals that 60,00 = vTTi(^ — e ~ Vi ) / P{ z )i with P(z) a quadratic 
polynomial in z. 

Identifying now the density with a new variable v/z we get 

- = p(z) = -5-r (u(z - it) - u(z + it)) ~ - artanh 6 ,oo = - log \ + ■ 
z 2m v ' z z 1 — B ,oo 

Exponentiating both sides we get an equation in two C* coordinates, z and x = e v . This 
defines a Riemann surface of genus one written in terms of two C* coordinates x and z, 
as expected for 5-dimensional Seiberg-Witten theory. The latter one can be brought to 
hyperelleptic form after the redefinition i— > x and a rescaling of the new variable x by 
P(z), which to the equation of the curve 



l 



x 2 



Yl(z-e^). (48) 

Three of the yiS are specified by the boundary conditions of the resolvent, and they depend 
also on the modulus t. To get the full solution, one more condition can be specified by 
imposing the A-period filling fraction. On the other hand, this A-period can be found from 
the knowledge of the above curve. 

Let us stress at this point, that the appearance of integrals (l47f ) of this form (in particu- 
lar leading to elliptic functions for SU(2) theory) is a consequence of the logarithmic branch 
cut structure of the derivative of the potential V(x) which appears in ( | 46 | ). Derivatives 
of both the quadratic and dilogarithm terms which appear in the potential are logarithms 
with branch cuts that we just need. This is therefore non-accidental that no other (special) 
functions arise in the potential. 

The above derivation could of course be generalized to arbitrary SU (n) theory, though 
technically the solution would be even more complicated. In this general case the well- 
known matrix model relations |30| for the filling fractions ai and the genus zero free energy 
F , written in terms of the rescaled variables 

If, dz dF Q I T dz 

a l = 7T- ( P logs — , -t— = - <t \ogX — , 

2m J A z dai 2 J B z 
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acquire a standard interpretation in Seiberg-Witten theory: they relate vevs of the Higgs 
field cti with the prepotential F , This completes then our solution of the Seiberg-Witten 
theory from the matrix model perspective. 

We also note that the matrix models which we consider are unstable due to the un- 
bounded potentials, as seen in figures Q] and El However this instability is not worse that 
the instability related to the condition of fixed filling fractions, which results in unequal 
level of a potential in each cut. Due to the condition of fixed filling fractions the eigenval- 
ues are prevented from falling into the infinite instability well in the same way as they are 
prevented from tunnelling between cuts. 

5 Relation to geometric engineering and topological 
string theory 

The purpose of this section is to give an unified geometric description of the 4- and 
5-dimensional supersymmetric gauge theories by geometric engineering of type II string 
theory and M-theory on non-compact Calabi-Yau spaces. This construction relates the 
matrix models discussed in this paper and the one in [221, [231 El]- For 4 dimensional 
supersymmetric gauge theories the subject started with the work of [T71 HHj and for the 5 
dimensional cases in [3 1 J . We give a short synopsis of these works focusing on the different 
geometries for the different Chern-Simons terms in the 5 dimensional case [31]. The main 
motivations are as follows. 

Physically the higher genus terms calculable in the matrix model are only relevant if the 
theory is embedded in gravity, which according to our present understanding is consistently 
possible only by geometric engineering of the corresponding string geometries. 

By moving in the moduli space, which is best described as the complex structure 
moduli space of the mirror geometry, one can connect the local matrix model descriptions 
of section [3] and [H to each other and to the Chern-Simons type matrix models that have 
been studied in [221 [231 [24] . 

The 5-dimensional theories are labelled by different integer Chern-Simons terms m G Z, 
which make an important difference in the identification with the matrix model, discussed 
in section 15.51 For this reason and because the absence of world-sheet instantons makes 
the theories simpler to describe, we discuss mainly theses cases [3 1 J . 

The non-compact geometries for the 5-dimensional case are the same as the ones used 
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in the 4-dimensional geometric engineering after a limit [18]. We can therefore focus 
in section [3] on the discussion of this limit, which leads to 4-dimensional spectral curve 
H(w, z\ Ha) = 0. We explain why it does not depend on m. 

For m = and ADE gauge theories G there is another point in the moduli space where 
one extracts the partition function of Chern-Simons theory on the lens space S 3 /Dg, where 
D G is a discrete subgroup of SU{2) acting on S 3 . The latter have a corresponding ADE 
classification and matrix model descriptions for these Chern-Simons theories have been 
suggested. The discussion in section EH is necessary to specify the general Chern-Simons 
matrix model limit in [23j. 

5.1 Geometric engineering of 5-dimensional gauge theories 

J\f = 1 supersymmetric gauge theory in 5 dimensions can be geometrically engineered 
from M-theory on certain non-compact Calabi-Yau threefolds M. The basic geometry 
of M is that of a resolved ADE singularity fibered over P 1 . It gives rise to a Af = 2 
supersymmetric gauge theory in five, and after a limit in four dimensions, whose gauge 
group G is the corresponding simply-laced group. The intersection of the rank(G) P 1 
curves in the fiber is the negative Cartan matrix of G. For non-simply-laced groups the 
fibration has to have non-trivial mondromy, see [31]. Matter can be added by additional 
blow ups [18j leading to singular fibres [31 J or by considering a higher genus base [HU [35]. 

The precise dictionary between manifolds M and the 5-dimensional gauge theories [HU 
l32l [33] contains a parameter m G Z, which labels a choice the Chern-Simons terms of 
the 5-dimensional gauge theory or equivalently the triple intersections of M. The choice 
m = corresponds to the matrix model discussed in section [H Since non-compact Calabi- 
Yau manifolds M on which M-theory compactifications give rise to 5-dimensional gauge 
theories with different Chern-Simons terms we need to describe how the differences in the 
non-compact mirror Calabi-Yau space W given by eq. ( I5"01) are reflected in the matrix 
model. At the end we propose matrix models describing topological string on a wide class 
of non-compact toric Calabi-Yau spaces, which correspond to the SU(n) cases. 

The different 5-dimensional theories, which however lead to the same gauge theory in 
4 dimensions, can be specified by a particular choice of the normalized triple intersection 
numbers of the non-compact Calabi-Yau, which appear in the 5-dimensional theory as 
Chern-Simons terms 




(49) 
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Here the latin indices label the vector multiplets. Note that (jl9i) is not gauge invariant, but 
it can be shown [36] that if M is spin with pi = 0, as it is the case if M is Calabi-Yau, the 
5-dimensional integral over (1491) shifts only by 27rZ. In M-theory on a Calabi-Yau M vector 
fields arise as A* = f e i.e. by reducing the three form field on a 2-cycle Q in a 
class in H 2 (M). Similarly the field strength arises as Fj = f c (dCp 4 \ so that (El9l) is just 
dimensional reduction of eleven dimensional Chern-Simons term f M A (dC)^ A (dC) ^ 
with Qjfc = J M A cjj A cjfe = Di fl Z?j fl .D^, where tUj are in H 2 (M) with compact support 
on d and £)j are dual divisors to [d]. Both definitions of the triple intersection need a 
regularization on non-compact Calabi-Yau spaces. The intersections are largely determined 
by the Cartan matrix of G and in [HI] it was shown that the discrete freedom in choosing 
different Chern-Simons terms is a discrete choice m G Z, which parametrizes that particular 
triple intersection mentioned above. For SU(n) all triple intersections will be made explicit 
in the prepotential below. 

5.2 Mirror symmetry for toric non-compact Calabi-Yau 

Mirror symmetry implies that the data of the 5 dimensional gauge theory as well as of 
the topological string on M are captured by a Riemann surface (j50p and a meromorphic 
differential A. One crucial point is that spectral curve of the 5-dimensional matrix model 
discussed in the last section is identified with that Riemann surface. More precisely in the 
moduli space of the Riemann surface the partition function of the 5-dimensional gauge 
theory emerges as a holomorphic limit of the topological string partition function near the 
large radius point \i a = 0. We review below how the mirror curve (1501 ) for the geometry 
described in (15.11 ) is constructed. 

The mirror manifold W of the manifold M is a conic bundle 

H(w, z; Ha) — uv = (50) 

branched over C* x C* at a complex family of mirror curves H(w, z; z a ) = 0, i.e. u, v G C, 
and w,z G C*. With the geometry comes a canonical meromorphic differential 

A = log(z)^. (51) 

From our perspective this differential will determine the filling fractions of the matrix 
model. 

For the SU(n) cases the well known construction of non-compact toric mirror symmetry 
applies [Hl[37], see [38] for a short review in the notation we use here. The manifolds M 
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and W are defined by charge vectors Qf G Z of the toric group action. More precisely 

M = {C k+3 - Z)/(C*) k , (52) 

where (C*) fe acts by Xi i— > Xj, a = 1, . . . , k on the complex coordinates Xj of C fe+3 with 
A Q G C*. The constraint Y^i=i Q? = 0) ^ a * s equivalent to C\(TM) =0. if is given by 
H = Yli=x x ii wherj^l xi G C* are homogeneous coordinates w.r.t. an additional C*-action 
and subject to the constraints 

fc+3 

(-lfoflxf = ^ aj Va. (53) 

-i=i 

Here fi ai a — 1, ... , h c °™ v , are complex moduli of (W, A), which are dual to the complexified 
Kahler parameters of M, while w and z in (l50ll denote the independent variables that 
remain after using the constraints (1531) as well as the additional C*-action on the Xi, 
mentioned above. H(w, z, // a ) = is central in the following discussion, because it is 
identified with the spectral curve of the matrix model. 

5.3 Geometric engineering of 5-dimensional SU(n) gauge theories 

Let us now give the dependence of the charge vectors Q1 defining M and W on m. 
The simplest case arises for SU (2) without matter. It was shown in [TH] that the possible 
5-dimensional geometry M correponds in this case to the anticanonical line bundle over 
the Hirzebruch surfaces F m = P(0 P i © C P i(m)), written here as the projectivization of 
line bundles over P 1 . The parameter m in this special case is precisely the m G Z freedom 
to choose different Chern-Simons terms. Recall that this non-compact Calabi-Yau spaces 
over the Hirzebruch surfaces F m are torically described by the 3-dimensional fan spanned 
by {ni, n 2 ,..., n 5 } = {(-1, 0, 1), (0, 0, 1)(1, 0, 1), (m, 1, 1), (1, -1, 1)}. The complete fan in 
the plane described by first two coordinates with the triangulation as shown in figure [4] 
defines torically the Hirzebruch surface F m . 

Each point rii corresponds to a divisor class £>i = {x^ = 0}, where are coordinates of 
the (l52l) . These divisors are not independent, but fullfill the relations Y^k=i n k,iDk for r = 5 
and i — 1, . . . ,d — 3. The triangulation specifies the intersection of M in the following 
way. DiDj — if % and j are not part of one triangle. DiDjDk = 1 if m,Tij,Tik span a 
triangle in the triangulation. Restricted to the Hirzebruch surfaces one has Y7k=i n k,iDk 

2 We use the same symbol Xi for the different coordinates of M and W. 
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F 



Figure 4: Toric fans for F, 



for r = 5 and i = 1, . . . , d = 2, D { fl Dj : = 1, i,j > 1 if rij, n,,- are on a lineal. There are two 
independent classes in F m , F = = D 5 and S = D 1 = S' — mF with 5" = D 3 . With 
the above information and the ring structure of the intersections it is easy to calculate all 
toric intersections. For the Hirzebruch surface FS = 1, F 2 = 0, S'S = and S 2 = —2. 
For M there is now a 2-dimensional complex dimensional compact divisor H = Di, which 
represents the Hirzebruch surface itself. H 3 = 4, HS 2 = —m, HF 2 = and HFS = 1, 
where S and S are non-compact divisors of M. Note that the triple intersections among 
the non-compact divisors corresponding to points on the boundary of the fan are not well 
defined. 

The charge vectors are coefficients of linear relations among the points YliiQi n i = 0- 
The one which do actually span the Mori cone are calculated using the secondary polytop 
to be Q 1 = (1, —2, 1, 0, 0) and Q 2 = (0, m — 2, —m, 1, 1). They correspond to curve classes 
and the intersection number of the a'th curve with the divisor D t is encoded in Qf. The 
Kahler cone is dual to the Mori cone and it is easy to see that it is spanned by F and 
H = S + (m + 1)F. 

These geometries are easily generalized to SU(n). To distinguish n odd and n even we 
define A(n) = |(1 — (—1)™). The points of the 3-dimensional fan and the Mori cone are 



3 We use hatted characters for non-compact divisors of M that restrict to compact divisors on F r , 
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Figure 5: Toric fans for the non-compact toric Calabi-Yau corresponding to 
SU (n) theory with the Chern-Simons term parameterized by m. 



One can see from the pictures that the inner points in the diagram correspond to toric 
compact divisors of the type in the non-compact threefold, which can be easily seen 
to be Hirzebruch surfaces. We denote these divisors H i7 i = 1, . . . , N — 1. From the left 
to the right Hi = F| m+2 j_ n |, % = 1, ... ,n — 1. With the description for the calculation of 
the intersection of the Hi it is easy to confirm [HI] that the cubic prepotential, which is 
not corrected by intantons for 5-dimensional theories is F — | \a% — %| 3 + m J2i a i — 

Sii,<a,t8 H h H i2 H i 3 n?=l( a v+l + • • • + O- 
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The standard parametrization of mirror symmetry using the Mori cone in (l53l) implies 
that /j, a = is a point of maximal unipotent monodromy and that near this point t a ~ 
2^ log(/i a ), with t a the compexified volumes of the curves corresponding to the Mori cone. 
For a — 1, . . . , n — 1 these curves are the fibres of the Hirzebruch surfaces. Geometrically 
their size is proportional to the mass of the W bosons in string units M string . It is given by 
the simple roots, i.e. t a = a a — a a+ i, a — 1, ... ,n— 1. The parameter t& :—t n corresponds 
to the complexified size of the base P 1 . The size appears as imaginary part while the 
S-field integral over the curve is the real part. 

5.4 4-dimensional Seiberg-Witten limit 

We come now to the engineering of the 4-dimensional gauge theory. The aim is to 
derive from the 5-dimensional spectral curve H(w, z; /i a ) = in a limit the 4-dimensional 
Seiberg-Witten (j^j) curve and show that A reduces to dS. Moreover in the 4-dimensional 
limit the differences given by the parameter m should vanish. While the 5-dimensional 
gauge theory discussed above is recovered in the large radius limit, the construction of 
the 4-dimensional gauge theory requires a double scaling limit in a region of the moduli 
space, called Seiberg-Witten point, where world-sheet instantons dominate. In fact it has 
been made explicit in [H] how the individual gauge theory instantons contributions of 
Seiberg-Witten gauge theory are encoded in a resummation of the world-sheet intantons. 

Non-abelian gauge groups in string theory have been first understood in the heterotic 
construction. For finding the Seiberg-Witten point it is instructive to recall the duality 
between the heterotic and type II strings. It states [17j that the heterotic dilaton S is 
geometrized and identified with j^t^, i.e. the size of the base of the fibration discussed in 
15.31 In the heterotic string one has Re(S') = 4? so that Im(t b ) = ^f. The asymptotic free 
g — > region of Seiberg-Witten must be at large volume of the base P 1 and hence near 
f'b '■= A*n = 0. However a double scaling limit, parametrized below by e — > must be taken 
to keep the mass of the W bosons finite, while taking the string scale to infinity. The 
mass of the W bosons is the only scale and the one-loop running of the gauge coupling in 
Af = 2 gauge theory implies \ ~ log (^f 1 )- Since VF-boson comes from branes wrapping 
the fibers one has U ~ C\og(tf), where tf refers to generic volumes of the fiber. The 
constant C can be determined from the dependence of the instanton contribution with 
instanton number k namely exp ^— ~ (jk^J ' wnere re f ers volumes of the fibers, 
i.e. masses of the W bosons, in gauge theory units. Putting the above facts together, the 
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double scaling limit is exp(t fo ) = (eA) 2n and tj = etf. In the B-model geometry, i.e. the 
manifold W, tf is mapped to period integrals over S , which come from periods over S 3 in 
the compact geometry. That is the Seiberg-Witten curve is near point where the divisor 
fi n = intersects n — 1 conifold divisors in the moduli space. The double scaling limit can 
also be viewed as resolving the point of tangency between the divisors of this intersection 
point, to obtain normal crossing divisors. 

Let us see for the SU{2) cases how (jjj) and dS is recovered in the limit from the 5d 
curve. According to the table in section [5731 and the description how to obtain H(w, z; fx a ) 
from (l53l we get 

u f z 2 ~ m 

H(w,z;jj a ) = w + ^- J -^ h 1 + z + fx b z 2 = (54) 

as the defining equation of the 5-dimensional curve. Here we identified in H = X^=i a; *> 
w := x 5 and z := x 2 . We set x 3 = 1 by the C*-action and eliminated x 4 and x\ by (153]) . Now 
we calculate the conifold discriminant and find that it behaves like A = (1 — 4fif) 2 + 0(fib), 
where all the m-dependence is in higher order in is proportional to A which and 

the scaling argument implies that the finite gauge theory coordinate u ~ B should appear 
as 

r/i = (1 -4///) = e 2 u, 

JJT b A 2 e~ s °/ 2 ( 55 ) 
V2 "(1-4 M/ )- u ' 

here e~ 5o//2 is an a priori undetermined constant and f l55ll are just the blow up variables 
to resolve the tangency between //& = and A = 0. An important point is that z is a 
C* variable in H(w, z; fi a ) = 0, but in flJJ) x is a C variable, which is clearly reflected by 
the log(z) dependence of A as opposed to the x dependence of dS. So we can think of f7]) 
as a compactification of H = in a special limit of its parameter space. To match the 
the differentials we must identify in leading order z = 1 — ex. Rescaling x = and 



eliminating the x term by x = x + ~ {^j= + ^y/wfj one gets independent of m in leading 
order w + = — x 2 — (l — jj^j . Rescaling further w — > w^/JIb and x — > l -^j^x and 



using (!55|) we get 

5.5 Matrix model description of the m^O cases 

For m 7^ the microscopic matrix model can be obtained starting with the modified 
Nekrasov partition function [33], which is written as in (fl~2l) with a linear modification and 
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quadratic modification in A by the second Casimir operator 



K\ = 



^Ai(Ai + l-2z) 



(56) 



and a linear term in A as 



■SU(n) 
5d 



2n\\\ e -mf3Y:M (l) \cn+hK x{l) z _ 



(57) 



A=(A( 1 ),...,A(")) 



If is straightforward to check [33], that this expression is reproduced by the applications 
of the topological vertex [39] to the geometries discussed in section 15.31 

Likewise it is simple to see, that if we apply the same steps as in section 0] to Z„ % 
we get fl39l) . except that we cannot view the whole expression, which follows from (15"?]) as 
a matrix integral, because of the different linear coupling (distinguished by ai) to n sets 
of eigenvalues in the potential. These linear terms are the only relevant change in the 
potential, since the shifts in the definition f|32l) give just rise to an overall constant and 
moreover the quadratic terms from the second Casimir operator will vanish in the H — > 
limit. 

In order to convert this to matrix model one has to introduce n matrices. The procedure 
is similar as in [23]. A linear term of the form ai Yli u f^ m a quadratic potential, which is 
integrated against the deformed measure can be absorbed into an interaction 



measure. This holds for the potential lf44|) . After viewing the n groups of u^ l \ I — 1, ... ,n 
as the eigenvalues of n matrices the interacting factors in the measure can be exponentiated 
to give a relatively complicated interaction term between the matrices [23]. The matrix 
models for m ^ describe, e.g. for Su(2) and m = 1 the anticanonical bundle over 
Fx, which by a blow down becomes P 2 . So it seems possible to write down a microscopic 
multimatrix model, which reproduces in a limit the simplest non-compact toric Calabi-Yau 



In this section we show that various well-known and seemingly unrelated matrix models 
arise in various limits of the 5-dimensional matrix model derived in section [H These 
interesting cases include Chern-Simons matrix models for lens spaces, as well as matrix 
models for line bundles over P 1 . Moreover, in the 4-dimensional limit the 5-dimensional 



between the different groups of eigenvalues 




e>(-3) -> p 2 . 



6 Relations to other matrix models 
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model reduces to the matrix model from section [3l which, as we discuss below, can be 
regarded as a generalization of the Eguchi-Yang matrix model (which itself corresponds 
just to the n — 1 case). In a sense, in this way we exhibit a part of a vast and unified 
matrix model landscape, whose various limits were introduced and considered in the past 
in an unrelated and independent fashion. We discuss these limiting cases in the following 
subsections. 

6.1 Chern-Simons matrix model for lens spaces 

The orbifold limit in the m = case is at the point in the moduli space where all tj, 
i = 1, . . . ,n, are small. In the complex structure variables this is near the vanishing of 
Wi — 1 — , i = 1, . . . , n — 1 and w n = , x r . Here is another matrix model 



known describing this topological string theory in the orbifold point [23]. Its spectral curve 

was analyzed in [24], [25] and claimed to be the curve H(x, y; z) = obtained from (l53l and 
the table in section 15.31 as described above. 

The matrix model presented in [23] looks almost the same as ours - the only difference 
is that its potential contains only the quadratic term. It does not contain neither the linear 
term in t, nor the dilogrithms. The absence of the linear term is easy to explain, as t = 
should in principle correspond to the orbifold limit. We also note that the dilogarithms 
are suppressed when their arguments have large negative values. It would be interesting 
to understand the relation between a continuation from the large radius to the orbifold 
point in more detail. We also note that slightly more general models for orbifold points 
corresponding to arbitrary lens spaces L(p,q) we recently found in (4TJ. They also have 
just the quadratic potential and the deformed measure, so belong to the same family of 
models as those described above. 

6.2 Matrix models for line bundles over P 1 

In [2T] topological strings for Calabi-Yau manifolds of the form 



were considered. Partition functions for such theories can be computed using the topolog- 
ical vertex [39] 





X p = 0(p - 2) © O(-p) -> P 




A 
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where the sum is performed over all two-dimensional partitions A, K\ = £\ Aj(Aj — 2i + l) is 
the second Casimir, q = e~ 9s , and t is the Kahler parameter. The vertex amplitude in this 
case Cr„ = Pq(X) 1 ' 2 is a square root of the deformed Plancherel measure (jSJ), so it leads 
to the same partition function as considered in [H]. The case p — 1 can be considered as 
the 5-dimensional U(l) theory, which is a special case n — 1 of our general result. It was 
shown in [21 J that the above partition function can be written as a one-cut matrix model 
with a deformed measure and the potential 

V{u) = (t-(p- l))u + \u 2 + 2 Li 2 (e"). 

For p = 1 this is indeed the potential which we obtain in n = 1 case. For arbitrary p only 
the linear and quadratic terms are modified (the latter due to the Casimir k\ in Zx p )- 

6.3 Eguchi-Yang matrix model 

We also note that in the limit of the 4-dimensional theory and for n = 1 our matrix 
model with the potential ( 1261 ) is closely related to the CP 1 matrix model introduced by 
Eguchi and Yang 

Z EY = J VMe -2NTrV^(M)^ 

OO OO * 1 

v ey = (_M + MlogM)-^^ liQ M J -^VM l (logM-^-). 

i=i i=i j=i ^ 



This matrix model was introduced in [191 ED], where it was also shown that its partition 
function Z EY reproduces the partition function of the topological CP 1 model considered in 
[40] , i.e. the topological sigma-model coupled to the two-dimensional gravity. This model 
has two observables denoted P = &o(P) and Q = o"o(Q), which correspond respectively 
to the identity and the Kahler class of P 1 . In addition, two infinite series of descendants 
cjj(P) and for % = 1,2, . . ., arise in this theory upon coupling to gravity. Coupling 

constants corresponding to these descendants are denoted by t^p and t^Q, and they appear 
in the Eguchi-Model as coefficients in two infinite series in the potential. 

We immediately see that a special case of the Eguchi-Yang model corresponding to 
ti,p = t it Q = for i > 1 is identical with the n = 1 case of our matrix model (l26l) . with 
our t identified with the coupling t ,Q m V EY . From physical point of view the vanishing 
of all coupling constants corresponds to replacing a matrix M by a chiral superfield [42] . 
From the model with n = 1 and vanishing coupling constants one can derive the limit 
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shape of partitions given by the arcsin-law [U, 43j . Our matrix model for Seiberg-Witten 
theory (l26l) is therefore a generalization of the Eguchi-Yang matrix model to arbitrary n. 
More generally, setting all t^p = but taking arbitrary t^Q is identical with the matrix 
model considered in [14], with arbitrary Casimirs taken into account. We note that some 
subtleties, adventages and drawbacks of the comparison of the sums over partitions arising 
from gauge theory instanton calculus, to the matrix model of the Eguchi-Yang type, are 
also discussed in [43] . 

7 Conclusions and discussion 

In this paper we provided a formulation of Seiberg-Witten theory in terms of matrix 
models. For 4-dimensional theories we derived 1-matrix models with a standard measure 
containing the Vandermonde determinant, while the 1-matrix models for 5-dimensional 
theories have a deformed measure. We also extended this approach to a matrix model 
formulation of topological string theory on Calabi-Yau manifolds which are related to 
Seiberg-Witten theories by geometric engineering. In general this involves multi-matrix 
models. 

The 1-matrix models which we propose were derived in two steps. Firstly, we introduced 
in the Nekrasov partition function an auxiliary parameter N, and rewrote this partition 
function as a matrix integral over matrices of size N. The matrix models obtained this 
way are rather non-standard 1-matrix models, with potentials which contain such special 
functions as logarithms of T functions. Then, the partition functions of the Seiberg-Witten 
theory were recovered in the 't Hooft limit. In this limit various simplifications take place 
and we end up with matrix potentials containing logarithmic and dilogarithmic terms. 

The derivation of these matrix models is just an initial step, which should be followed 
by a detailed analysis of their properties. We use the Migdal-Muskhelishvili formula for 
the resolvent to check that the 1-matrix models for the 5-dimensional theories yield the 
expected spectral curves with the correct differential A defining the filling fractions. The 
4-dimensional spectral curves and the corresponding differential dS are reproduced by in 
the gauge theory limit of geometric engineering. It could have also been derived with the 
Migdal-Muskhelishvili formalism for the resolvent in 4-dimensional theory directly. 

The general framework for solving matrix models based just on the underlying spectral 
curve and the differential, which was developed by Eynard and Orantin [HI 01] and checked 
in very similar situations in [45l [46] , should give the higher genus, as well as open amplitudes 
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in all our cases The holomorphic anomaly equations and the gap condition applies 
likewise to the 4-dimensional cases [12] and to the 5-dimensional models [38j. It provides 
the most efficient way to calculate the closed higher genus amplitudes. 

Certainly a perturbative expansion of the matrix models discussed here is possible. Of 
course we hope that the knowledge of the microscopic matrix model action will provide a 
route towards better understanding non-perturbative effects in Seiberg-Witten theory and 
topological strings. Techniques of analyzing non-perturbative effects in matrix models, 
being currently developed by [47l [48], require the knowledge of the potential, which we 
provide in this paper. They would also be of great interest from the point of view of black 
hole physics and OSV conjecture [19], and complementary to the proposition discussed in 
PI EH- 

An important aspect of the matrix models derived here is their relation to other families 
of matrix models. Firstly, matrix models for the SU(n) 4-dimensional Seiberg-Witten 
theory can be regarded as a generalization of the the Eguchi-Yang matrix model (which 
corresponds to n = 1). Higher Casimir operators, which we did not analyse in this paper, 
but in the context of n = 1 model they were taken into account into [HI [43], would 
correspond to gravitational descendants in the Eguchi-Yang model. In the same spirit, 
matrix models for 5-dimensional SU(n) theory derived here are generalizations of the 
matrix model for bundles over P 1 considered in [2T] (which also corresponds to n — 1). It 
is also reassuring that the models which we derived contain no other terms than logarithms 
and dilogarithms, which also appeared in [21]. While our models, as well as those in [2"T] , 
correspond to the large radius point in the moduli space, it is also possible to consider 
other regions of the moduli space. In particular, in the limit where dilogarithms in the 
potential for the 5-dimensional theory are suppressed, we recovered Chern-Simons matrix 
models for lens spaces which correspond to orbifold points in the moduli space [23 l [241 [25] . 
The similarity of the microscopic description at the large radius point and at the orbifold 
point (just the suppression of dilogarithmic terms) is surprising in view of highly involved 
transformations of partition functions between these points. 

The above results indicate a possible unified microscopic description of the multitude 
of matrix models, which describe various regions in the moduli space. Similar ideas of the 
unification of matrix models to some overall matrix-model- M-theory (in analogy to string 
theory) were advocated in [52, 53J. 

Our results should also provide a direct link between fermionic formulations of matrix 



4 Some checks in 4-dimensional Seiberg-Witten theory have been made in |13j . 
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models, gauge theories and topological strings. The fermionic viewpoint on matrix models, 
intimately related to conformal field theory techniques, was originally advanced by Kostov 
|54j . Recently the (matrix model) formalism of Eynard and Orantin was also rephrased in 
a fermionic language [55]. On the other hand, topological strings and gauge theories can 
be formulated in terms of two-dimensional chiral fermions too [TU1 [39l EU E7] EH [59]. Such 
fermions live on curves which, in the present context, are simultaneously matrix model 
spectral curves and Seiberg-Witten curves. This convinces that all those fermions have a 
common origin. Moreover, in [57] a formalism was presented which relates both matrix 
models and Seiberg-Witten theory to P-modules. Finding how matrix models presented 
in this paper fit into this formalism would certainly be of interest. 

Let us also stress that the idea of expressing gauge theory quantities in terms of matrix 
integrals has been explored to much extent in the context of Af = 1 theories within the 
so-called Dijkgraaf-Vafa formalism [60] ETj . Af = 1 theories can be treated as deformations 
of Af = 2 theories, and it was also shown that certain scaling of Dijkgraaf-Vafa matrix 
models leads to the Af = 2 answers (HJ [62]. The difference with the case considered in 
this paper is that now we get the answer which gives Af = 2 results without a need for any 
rescaling. Some other matrix models or related ideas corresponding to gauge theories were 
also proposed in [63], EH EHl [66], ETj. It would be interesting to relate these points of view. 
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Appendix - elliptic integrals 



To derive the spectral curve for 5-dimensional SU(2) theory we need to perform the 
following integral 

v^ 1 (z — v)~ x dv 



y/(v — a) (v — b) (v — c) (v — d) 



9l F(y(v),k) + g 2 U(y(v), 5(0), k) + g 3 Tl(y(v), 5(z), k), 



which is expressed as a combination of the elliptic functions, with the following notation. 
The coefficients in this combination are equal 

2 2(6 -a) 

9i = 1 92 = , = , (58) 

a(a — z) y/(b — c) (a — a) zab^J (b — c) (a — a) 

2(6 -a) 

9s - 



z(z — a)(z — b)y/(b — c)(a — d) 
while the arguments of F and II are given as 

(a-d){b-v) 6(z) _J_ 2 1 (b-d)(a-c) 

V[V) ~ (b - d)(a - vY [) ~y(zY ~ y(c)~ (a-d)(b-cY [ ™ } 

The elliptic functions F and II, of the first and the third kind respectively, are 

,n f v dt P ,x F v dt 

F(y,k)= __=__=_, U(y,5,k) 



lo v /(l-t2)(l_ fc 2 t 2)' ^' ' 7 (l-St 2 )y/(l-t 2 )(l-kH 2 )' 

(60) 

The complete elliptic integrals are defined as K(k) = F(l, k) and 11(5, k) = 11(1, 5, k). Of 
particular importance in our computation are the following addition formulas [68] 



F(sin$,£;) ±F(sin/3,fc) = F(sin^,A;), (61) 
n(sini?,5,fc) ±n(sin^,5,fc) = Il(siny?, 5,k) ± J— — — — artanhG, (62) 



(5-l)(5-k 2 ) 



where 



@ _ sin^ sin/3 siny? -^(fl - l)(5 - k 2 ) 

1 — 5 sin 2 if + 5 sin $ sin /3 cos cp yl — k 2 sin 2 (p 



and the angle ip is determined as 



cos "d cos /3 =F sin -d sin /3 >/(l — k 2 sin 2 $) (1 — k 2 sin 2 /3) 
1 — k 2 sin 2 1? sin 2 /3 
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